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Abstract 

Recent innovations on the differential calculus for functions of non- 
commuting variables, begun for a quaternionic variable, are now ex- 
tended to the case of a general matrix over the complex numbers. The 
expansion of F(X+Delta) is given to first order in Delta for general 
matrix variables X and Delta that do not commute with each other. 
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1 Introduction 

In a recent paper [1] I showed how to expand 

F(x + 8) = F(x) + F (1) (x) + 0{8 2 ) (1.1) 

when both x and 8 were general quaternionic variables, thus did not commute 
with each other: 

F^ 1 \x) = F'(x)5 1 + [F(x)-F(x*)}(x-x*)- 1 ^ 5 = 6 l + 8 2 , (1.2) 

with specific formulas on how to construct the two components of 8. 

Now we shall extend that analysis to a more general situation. 

Consider the N x N matrices X over the complex numbers and arbitrary 
analytic functions F(X) with such a matrix as its variable. We seek a general 
construction for the first order term F^(X) when we expand F(X+A) given 
that A is small but still a general N x N matrix that does not commute with 
X. 

The first step, as before, is to represent the function F as a Fourier 
transform, 

F(X) = J dpf(p)e^ x (1.3) 

where the integral may go along any specified contour in the complex p-plane; 
and then we also make use of the expansion, 

e (X+A) =e X [1+ f 1 dse sX A e sX + 0(A 2^ (L4) 
JO 

Another well-known expansion, relevant to what we see in (11.41) . is 

e A B e~ A = B + [A, B] + [A, B]} + I [A, [A, [A, B]]} + ..., (1.5) 
involving repeated use of the commutators, [A, B] = AB — BA. 

2 Diagonalization 

The first step is to assume that we can find a matrix S that will diagonalize 
the matrix X at any given point in the space of such matrices. 

A = SXS-\ A^ = 8 %J \i, i,j = l,...,N (2.1) 
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and we carry out the same transformation on the matrix A: 

B = S A S- 1 



(2.2) 



but, of course, the matrix B will not be diagonal. 

Our task is to separate the matrix B into separate parts, each of which 
will behave simply in the expansion of Eq. (11.51) . The first step is to recognize 
that the diagonal part of the matrix B, call it B , commutes with the diagonal 
matrix A and thus we have 

e tA B e~ tA = B , (2.3) 

where we have t = —isp from (11.31) . (11.41) . 

Next we look separately at each off-diagonal element of the matrix B: that 
is, Buj\ is the matrix that has only that one off-diagonal (i ^ j) element 
of the given matrix B, and all the rest are zeros. The first commutator is 
simply 

[A, B (i>j) ) = m B (i>j) , Ti j = Xi - Xj (2.4) 
and then the whole series can be summed: 

e tA B {hj) e- tA = e tr " B {iJ) . (2.5) 

Putting this all together, we have 

S S- 1 = f dpfip) ip e^ A f 1 ds [B Q + J2 e- isWl3 %,-)]• (2-6) 

It is trivial to carry out the integrals over s; and we thus come to the 
final answer 

F«(X) = F\X) A + Y\F{X) - F{X - nj I)\ A (m) (2.7) 

where 

A M = 5- 1 B^ S, n = 0,(i,j). (2.8) 

3 Discussion 

The general structure of the result, Eq. fl2.7[) . is similar to what we found 
in earlier work, Eq. (11.21) : the first term F'(X) looks like ordinary differen- 
tial calculus and goes with that part of A that commutes with the local 
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coordinate X; the remaining terms are non-local, involving the function F 
evaluated at discrete points separated from X by specific multiples of the 
unit matrix I. 

While this final formula appears not as explicit as the previous result 
found for quaternionic variables (or for variables based upon the algebra 
of SU(2)), in any practical situation we have computer programs that can 
calculate the matrix operations referred to above with great efficiency 

The quantities r^, which may be real or complex numbers, can be called 
the "roots" following the familiar treatment of Lie algebras. They have some 
properties, such as = — rji and sum rules that involve traces of the original 
matrix X and powers of X. 

What happens if the eigenvalues of X are degenerate? Suppose, for ex- 
ample, that Ai = A2. This means that ri2 and r%\ are zero. If we look at Eq. 
(12.71) . we see that the terms A(i >2 ) and A( 2) i) then have the coefficient F'(x). 
Thus they simply add in with Ao- In the extreme case when all the eigenval- 
ues are identical, then the answer is F(X + A) = F(X) + F'(X) A + 0(A 2 ), 
which is old fashioned differential calculus for commuting variables. 

Appendix A 

A more abstract form of the result Eq. (12. 7p is the following. 



Here the matrix C is defined implicitly, through its commutator with X, 
rather than explicitly; and the matrix A is the same as previously discussed. 
One may readily confirm the correctness of this formula in the case of F(X) = 



This alternative formalism may also be applied to the case of a quater- 
nionic variable x, which was studied in reference [I]. In that case we find 
A = S x and C = pr[x,5]. 

It is interesting that this new formalism manages to hide the non-locality, 
which was a prominent feature of the original analysis. 



F^(X) = F'(X) A +[C,F(X)\, 
[C,X} = A-A , [A ,A] = 0. 



(A.l) 
(A.2) 



X n . 
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